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INPUT: Matrix P 
Constant c 
Parameter k 
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LET M = (I -(l-c)P)/c 



I 
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COMPUTE d = |A/| 2 



508 



COMPUTE U k = k eigenvectors 
of (dl- MM 1 ) by calling 
kEigen(M, n, k, d) 
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COMPUTE V k = ^eigenvectors 
(dl - M T M) by calling 
kEigenfM 7 , n, k, d) 
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LET matrix D contain eigenvalues 
of (dl - MM 7 ); let S = (dJ- D) 1/2 
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RETURN (U„ S, VJ 
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INPUT: 
Matrix M, Dimension dim 
Parameters k, d 
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CHOOSE k vectors a p . . ., a k ^ 
randomly from R dim ; form matrix 
A with columns a p .... a k , 
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COMPUTE B, with column B t 
computed by Mult(M, a f d) 
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SET .4 = Orthonorm(B) 
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Mult 



INPUT: Matrix R 
Vector v 
Constant d 
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COMPUTE x = R T v 
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COMPUTE y = Rx 
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COMPUTE & RETURN (dv-y) 



